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Abstract. By analytically solving the master equation, we investigate quantum state transfer, creation and 
distribution of entanglement in the model of Milburn's intrinsic decoherence. Our results reveal that the 
ideal spin channels will be destroyed by the intrinsic decoherence environment, and the detrimental effects 
become severe as the decoherence rate 7 and the spin chain length A*' increase. For infinite evolution time, 
both the state transfer fidelity and the concurrence of the created and distributed entanglement approach 
steady state values, which are independent of the decoherence rate 7 and decrease as the spin chain length 
A'^ increases. Finally, we present two modified spin chains which may serve as near perfect spin channels 
for long distance state transfer even in the presence of intrinsic decoherence environments T[p{t)]. 

PACS. 03.67.-a Quantum information - 03.67.Mn Entanglement production, characterization, and ma- 
nipulation - 03.65.Yz Decoherence; open systems; quantum statistical methods 



1 introduction 

In quantum information processing (QIP), it is desirable 
to achieve a high-fidelity transfer of quantum states be- 
tween different parts, such as the core processor, stor- 
age, etc., of a quantum computer. To this end, a vari- 
ety of solid-state spin networks with always-on interac- 
tions have been proposed [1-16]. Particularly, Christandl 
et al. showed that with elaborately designed modulated 
exchange couplings between neighboring spins, one can 
implement perfect quantum state transfer (QST) over ar- 
bitrary distances between the opposite ends of a XX spin 
chain or between the two antipodes of the one- link and the 
two-link hypercubes with however the maximum perfect 
communication distance 21og3A [3,4]. In addition, they 
also showed that these modulated spin structures can dis- 
tribute arbitrary entanglement between two distant par- 
ties. Zhang and Long et al. [13] realized this perfect state 
transfer algorithm in a three-qubit XX chain using liquid 
NMR system. Later, Shi et al. presented a class of more 
general pre-engineered perfect spin channels [6] according 
to the spectrum-parity-matching condition (SPMC) they 
deduced. Then, Kostak et al. [14] established a general 
formalism for engineering spin Hamiltonians for perfect 
state transfer in networks of arbitrary topology and cou- 
pling configuration. Christandl's innovative works were ex- 
tended by Jafarizadeh and Sufiani in a recent work [15], in 
which they adopted distance-regular graphs as spin net- 
works and found that any such network (not just the hy- 
percube) can achieve unit fidelity of state transfer over 
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arbitrarily long distances. Moreover, D'Amico et al. [10] 
showed that one can create and distribute entanglement 
with an interaction-modulated Y-shaped spin network, 
particularly, with a slightly complicated bifurcation struc- 
ture, the distributed entanglement can be frozen when a 
phase flip is applied to one spin out of each pair. 

In addition to the above-mentioned protocols which 
mainly concentrated on spin chains with nearest-neighbor 
(NN) couplings, in Ref. [17] Paternostro et al. studied QST 
in imperfect artificial spin networks with all the qubits 
are mutually coupled (in which the usually assumed NN 
coupling is invalid). They presented a strategy to avoid 
the spoiling effects of these redundant connections with 
a modification of the couplings of the first and the last 
qubits in the chain, which enables nearly optimal state 
transfer. Then in Ref. [18] Kay demonstrated that per- 
fect state transfer is also possible in the presence of next- 
nearest-neighbor (NNN) couplings. Moreover, compared 
to the case where the system contains only two-spin inter- 
actions, the authors in Ref. [19] presented a scheme of QST 
by introducing the three-spin interaction, and showed that 
they can significantly increase the speed of QST in an 
XY chain. Besides the spin-half systems, state and en- 
tanglement transfer driven by a bilinear-biquadratic (BB) 
spin-1 Heisenberg chain was also discussed recently [20], in 
which the authors concentrated on the relations between 
the transfer efficiency and the quantum phase transitions. 

Most recently, a milestone work appears in Ref. [21] 
presented a control-limited scheme [22] for perfect state 
transfer through a pre-engineered spin chain with the help 
of local end-chain single-qubit operations. While nearly 
all of the previous schemes whose achievements of perfect 
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state transfer relies crucially on the preparation of the spin 
medium in a fiducial pure state, the authors in Ref. [21] 
demonstrated that state initialization of the spin medium 
is inessential to the performance of the protocol if proper 
encoding at the end of the chain is performed. The key re- 
quirements for their scheme are the arrangement of proper 
time evolution and the performance of clean projective 
measurements on the two end spins. This innovative work 
considerably relaxes the prerequisites for obtaining reli- 
able QST across interacting-spin systems. Stimulated by 
this innovative work, in Ref. [23] Markiewicz and Wiesniak 
proposed a special type of two-qubit encoding strategy for 
perfect state transfer, where no remote-cooperated global 
state initialization and any additional communication are 
needed. 

Apart from these exciting progresses, we noted that al- 
though there are several works [24-36] concerning the de- 
coherence effects on entanglement dynamics, stiidics thus 
far has seldom consider the influence of different kinds of 
decoherence scenarios on transfer of quantum states due 
to the complex and unclear mechanism of its interaction 
with the environments. However, from a practical point of 
view, all the real physical systems, especially a solid-state 
system, will unavoidably be influenced by its surrounding 
environments. This influence can cause the initial state 
of the system of interest becomes entangled with the en- 
vironment in an uncontrollable way, and it is just this 
entanglement of the system with the environment that 
causes decoherence. The decoherence can greatly affects 
the transfer efflciency of quantum states, as well as gener- 
ation and distribution of entanglement, and thus becomes 
one of the dominating obstacles baffling the physical im- 
plementation of QIP. It is therefore of great importance 
and fundamentally interesting to find ways to prevent or 
minimize the detrimental effects in the practical realiza- 
tion of QIR 

The standard way to investigate decoherence is to con- 
sider the system of interest as a part of a larger closed 
system involving the environment, and the density oper- 
ator of the system can then be obtained by tracing out 
all other degrees except quantum states of the system. In 
the present paper, however, we would like to resort to a 
different approach, i.e., the scenario of the so-called intrin- 
sic decoherence proposed by Milburn [37], who modified 
the Schrodinger equation in such a way that quantum co- 
herence is automatically destroyed as the system evolves. 
Such a consideration is fed by two motivations. First, this 
model is amenable to exact analytical treatment as we will 
see, one can determine the density operator of the system 
at arbitrary time t by the sole knowledge of the eigenvalues 
and eigenvectors of the system. Second, although the ab- 
sence of unitarity for a closed system in this model makes 
it unlikely to be a fundamental description of decoherence, 
its stochastic behavior in time evolution may still be an 
effective approximation for describing the phenomenon of 
the system. For example, it has been applied to describe 
decoherence of a single trapped ion due to intensity and 
phase fluctuations in the exciting laser pulses [38] . Dynam- 
ics of the mutual entropy of two-coupled Josephson charge 



qubits with intrinsic decohercincc; has also been studied re- 
cently [39]. Moreover, as pointed by the authors of Ref. 
[29,40], this model may be available in approximately de- 
scribing the non-dissipative decoherence of several physi- 
cal systems in the presence of white noise. 

2 General formalism 

In this paper, we consider quantum state transfer proper- 
ties in the model of Milburn's intrinsic decoherence [37]. 
The kernel of this decoherence scenario is the postulate 
that on sufficiently short time steps the system does not 
evolve continuously under unitary evolution but rather in 
a stochastic sequence of identical unitary transformation, 
which can account for the disappearance of quantum co- 
herence as the system evolves. Based on this assumption, 
Milburn obtained the master equation (in units of H) gov- 
erning the time evolution of the system 

^ = - [exp{-i-fH)pexp{ijH) - p], (1) 
at 7 

where 7 is the intrinsic decoherence parameter (the mean 
unitary time step). Expanding Eq. (1) to the first order in 
7, one finds 

^ = -^[H,p]-l[H,[H,p]]. (2) 

The first term on the right-hand side of Eq. (2) gen- 
erates a coherent unitary time evolution of the system, 
while the second term, which does not commiite with the 
Hamiltonian, represents the decoherence effect on the sys- 
tem and generates an incoherent dynamics of the system. 
In the limit of 7 ^ 0, the ordinary Schrodinger equation 
is recovered. 

To solve Eq. (2), one can define three auxiliary super- 
operators J, S and L, which satisfy 

Jp = lHpH, Sp = -i[H,p], Lp = -^{H\p]. (3) 
Prom Eq. (3) it is straightforward to show that 

exp(Jr)p(t) = )^H'p{t)H\ 

exp{ST)p{t) = cxp{—iHT)p{t) cxp{iHT), 
exp{LT)p{t} = exp (- yi?') p{t)exp {'^H^) ■ (4) 

Thus Eq. (2) simplifies to dp/dt = {J -\- S + L)p, and 
its formal solution can be written in terms of the Kraus 
operators Mi{t) as 

oo 

p{t) = j2m)p{mi{t), (5) 

1=0 

where p{0) denotes the initial state of the system, Mi{t) = 

(jty/'^H'^ cxp{-iHt) cxp(-7ii?2/2)/VZ! satisfies the rela- 
tion Y,Zo ^i{t)Mi{t) = 1 for all time t. 
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If we rewrite p{0) in forms of the energy eigenstate 
basis as p(0) = J2kk' '^kk'\''Pk}{''Pk'\, then we obtain 



Pit) 



E 

kk' 



Qkk' ex.p[-it{Ek-Ek')- — {Ek-Ek'Y] \'^k){'4'k' \ , 



(6) 

where akk' = {'>Pk\p{0)\ipk') , Ek and \tpk) are eigenvalue 
and the corresponding eigenvector of the considered sys- 
tem. 

For the special case that p(0) is an eigenstate of the 
system, Ukk' only when k = k', Thus from Eq. (6) 
one can obtain p{t) = akkli'k) {tpk\ = p(0)i the system 
will be unaffected by the intrinsic decoherence during the 
time evolution process. 

Furthermore, for a spin chain Hamiltonian commutes 
with the total z component of the spin, i.e., [H,a^^^] = 
0, where a^^^ = J^i'^i^ the 2^ (g) 2^ Hilbert space can 
be decomposed into N + 1 different invariant subspaces, 
each of which is a distinct eigenspace of the operator a^^^, 
and a system prepared in these subspaces will remains 
in them. In the single-excitation invariant siibspace 'Hi 
spanned by the site basis |n) = cr+|0)®^ (n = 1, 2, N), 
one can rewrite p{t) as p{t) = ^nm^nm\'n){m\, then in 
the standard basis {|00), |01), |10), the single qubit 

reduced density matrix can be obtained as 



Pi{t) 



-hi 

bu 



(7) 



Similarly, one can obtain the two-qubit reduced den- 
sity matrix between qubits i and j as 



1 



Pij{t) 



Hr-bjj 0^ 
bjj bj, 
bij bu 



(8) 







0, 



In this paper, we use the fidelity F = {tjj{0)\pi{t)\ip{0)) 
as an estimation of the quality of the state transfer from 
the sender to the destination qubits [1], and adopt the 
concept of concurrence C = max{0,Ai — A2 — A3 — A4} 
as a measure of the pairwise entanglement [31,32]. Here 
the quantities Ai {i = 1, 2, 3, 4) are the square roots of the 
eigenvalues of the product matrix R = p{ay (g) (T^)/9*(cr^ (g 
cr^) in decreasing order. 

From Eqs. (7), (8) and the above definitions about 
transfer fidelity and concurrence, one can obtain directly 
that F{N,t) = \bNN\ and Cij{N,t) = 2|6y| for a state 
initially prepared in the A^-dimensional subspace Hi- 

Another quantity related to the efficiency of the quan- 
tum spin channel of interest is the fidelity averaged over 
all pure states in the Bloch sphere. The state of the whole 
system at the initial time t = can be written as 



|V(0))=cos-|0)+e'*sin-|s), 



(9) 



where |0) = |00...0), \s) = a+\0)®^ , 6 and (j) are arbitrary 
phase angles. 



For this type of initial state, its dynamics is completely 
determined by the evolution in the zero and single exci- 
tation subspace Hqqi. From Eq. (6) one can obtain the 
state at time t as 



N 



pit) = cos^ ^\0){0\ + sin^ - ^ anm\n){m\ 

n,m=l 

( ^ \ 

+ sin - cos - bn\n){Q\ + H.c. j 
with the coefficients a„m and 6„ given by 



(10) 



N 



^ ^ ^k,s^k' ,s^k,n^k' ,m ^ 

k,k' = l 



exp 

TV 



-it{Ek — Ek') — ^(-^fe ~ Ek')^ 



( it \ 

bn = ^Cfe,sCfe,„ exp i^itEk - y-^^ftj > 



(11) 



where Ck,n is the amplitude of coefficient for the state |n) 

in the eigenstate \k) = Yl,n=i Cfe,n|")- 

Then by tracing off the states of all other spins except 
i from p{t), one has 

«w = (A;.""?"j/'':*ilr^). (12) 



6,6**' sin f cos \ 



an sm 



Prom Eqs. (9), (12), the fidelity F = (V'(0)|pi(t)|V'(0)) 
can be obtained as 



F = cos^ -\ \ — an aiv? - + 2|6j| sin^ - cos a 



an sm 



(13) 



where a = arg(6j) denotes the argument of the complex 

number bi. 

Thus the average fidelity F = {tf;{0)\pi{t)\ilj{0))df2 
can be calculated as 



=, |6i|cos(a) an I 

r = -r- r ~. 

3 6 2 



(14) 



From Eq. (11) one can see that in the absence of in- 
trinsic decoherence (i.e., 7 = 0), the equality an = 
holds, thus Eq. (14) reduces to Eq. (6) in Ref. [1], which 
describes average fidelity in the non-disturbed case. 



3 State transfer in decoherence spin channels 

We first consider quantum state transfer via spin chain 
governed by the XX Hamiltonian 



J 



N-l 



H=-J2^«+l+''>n+l), 



(15) 
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where a" (a = y, z) are the usual Pauh matrices of the 
nth qubit. 

For this model, Christandl et al. have shown that per- 
fect state transfer from one end of the chain to another 
is only possible for the case of chain length N = 2 and 

= 3, respectively [3,4]. Here we show that this ideal 
communication channel will be destroyed under the influ- 
ence of intrinsic decoherence. 

The eigenvalues and eigenvectors of the Hamiltonian 
(15) can be obtained as 



El. — 2 J cos 



nk 
N + 1' 

— N 



1^' 



nkr 



N+1 



(16) 



We first consider transfer of an excitation across the 
chain. For this purpose, we assume the system is initially 
prepared in the state |no). In the energy eigenstate basis, 
|no) can be expressed as 



Thus one has 



N 



N- 



sm 



k=l 



Trkno 
N+1 



\k). 



(17) 



N 



k,k' — l 



(18) 



Combination of Eqs. (6) and (18) gives rise to 



N 



N 



Pit) 



E E 

n,7n— 1 k,k' — l 



T:kn \ , ( T:k'm 
sm I sm 



N + 1 



N + 1 



nknQ \ ( ■Kk'riQ 
X sm I — - — - sm — - — - I X 



exp 



exp 



N +1 



-i2Jt cos 



N +1 



-2J^-ft cos 



■nk 
N +1 

irk 



cos ■ 



Trfc' 
N + 1 

nk' 



N + 1 



N+1 



\n) (m|. 
(19) 



For initial state |1) prepared in the input node A, the 
transfer fidelity of the output state in node B can be ob- 
tained from Eq. (19), and typical plots for the cases of 
N = 2 and TV = 3 with different decoherence rates are 
shown in Fig. 1, where the coupling constant J is chosen 
to be 1. In big contrast to the ideal case (i.e., 7 = 0), one 
can see that the transfer fidelity F behaves as a damped 
oscillation as the time t evolves. This phenomenon can be 
understood from Eq. (19), where the product of the first 
five terms on the right-hand side causes the oscillations, 
and the last term introduces the amplitude damping. With 
the increase of the decoherence rate 7, or the chain length 
N, the detrimental effects becomes more severe and there- 
fore more quantum state information will be lost. Thus for 
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Fig. 1. (Color online) Dynamics of the state transfer fidelity F 
for the XX chain with identical interactions. The decoherence 
rate is given by 7 = 0.1 (black), 7 = 0.2 (red) and 7 = 0.3 
(blue). 



spin networks with identical neighboring qubit couplings, 
even if for the one-link and two- link hypercube geometries, 
perfect transfer of an excitation is still impossible in the 
intrinsic decoherence environments. 

For infinite time t, the system evolves into a steady 
state with the transfer fidelity arrives at an asymptotic 
value F'**®'*'^^(iV), which can be obtained by combination 
of Eqs. (7), (19) and taking the infinite-time limit. After 
a tedious computation, we obtain 



2(iV+ 1)' 



(20) 



Clearly, this steady state transfer fidelity is indepen- 
dent of the decoherence rate 7, and it solely decreases with 
the increase of the chain length N. 

Next we consider time-dependence of the average fi- 
delity for the XX spin chain with identical interactions 
and subject to intrinsic decoherence environments, with 
initial state prepared in the form of Eq. (9) in node A, 
i.e., s = 1. From Fig. 2 one can see clearly that the av- 
erage fidelity F also behaves as a damped oscillation as 
the time evolves. Here the relative small value for the case 
of iV = 2 is due to the fact that the phase of the state 
at node B is uncorrected, i.e., a is not a multiple of 27r. 
When f — )• 00, the average fidelity also arrives at a steady 
state value, which is independent of the decoherence rate 
7, and can be obtained analytically by taking the infinite- 
time limit of F from Eqs. (11), (14), and (16) as 



^steady i^^-j = 



6iV+ 17 

12iV+ 12 
2iV + 3 



4A^ + 4 



if N e odd, 



if N ^ even. 



(21) 



Contrary to that of the initial state |1), this steady 
value does not decrease monotonously with the increase 
of the chain length N. However, as can be seen from Eq. 
(21), they decrease with the increase of the odd and even 
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Fig. 2. (Color online) Dynamics of the average fidelity F for 
the XX spin chain with identical interactions, where the deco- 
herence rate is given by 7 = 0.1. 



A'', respectively, and approach to the asymptotic value 0.5 
in the limit oi N ^ 00. 

In the following we discuss quantum state transfer in 
intrinsic decoherence spin channels with fixed but differ- 
ent couplings between qubits. We consider the following 
modified Hamiltonian 



N-l 

^ = 2^ 7) \^n<^n+\ 



1), 



(22) 



where Jn,n+i = A\/ n{N — n) is the modulated exchange 
coupling, and A is a scaling constant. 

The above Hamiltonian is identical to the represen- 
tation of the Hamiltonian Hs of a fictitious spin S = 
(N — l)/2 particle: Hs — AS'^, where Sx is its angu- 
lar momentum operator in x-direction and A is a scaling 
constant. For this Hamiltonian, its eigenvalues and corre- 
sponding eigenvectors can be obtained as [43] 



N 



Ek = i-N + 2k^l)X, |fc) = ^cfc,„| 



(23) 



n=l 



where the coefficient Ck.n is given by the following recur- 
sion relations 



ci,i = l/2(^-i)/2 



^[n-l){N -n + l) 



{n ^ 2). 
(24) 



For this modulated chain, it has been shown that one 
can achieve perfect state transfer between the input node 
n and the output node N —n-\-l after a time = 7r/2A and 
at intervals of tt/A thereafter in the absence of decoher- 
ence environment [3,4]. When the intrinsic decoherence is 
present, however, this ideal spin channel will be destroyed, 
and it acts as an amplitude damping quantum channel as 
the rescaled time \t evolves. As can be seen from Fig. 
3, the transfer fidelity F oscillates around a steady state 
value, with the amplitude decreases gradually. This detri- 
mental effects becomes more and more severe with the in- 
crease of the decoherence rate and the spin chain length. 




Fig. 3. (Color online) Dynamics of the state transfer fidelity F 
for the XX chain with modulated interactions, (a) chain length 
A'^ = 100 with different decoherence rate; (b) decoherence rate 
7 = 0.15 with different chain length. 



which is in consistent with the cases of the two- and three- 
site spin chains with identical interactions (In fact, they 
are two special cases of the interaction-modulated spin 
chain). This puts new constraints on these spin chains for 
long distance quantum state transfer. When t — > 00, the 
transfer fidelity reaches a steady state value, which can be 
obtained from Eqs. (6), (7), (23), and (24) as 



F 



(Steady 



{N) = 



1 



22N- 



N 

n 

k=2 



4- 



1 



(25) 



The steady state transfer fidelity of the interaction- 
modulated spin chain is still independent of the decoher- 
ence rate 7, and its magnitude is larger than its unmodu- 
lated counterparts [cf. Eqs. (20) and (25)], thought it still 
decreases with the increase of the chain length N . 

On the other hand, since the detrimental effects be- 
come severe as the rescaled time \t evolves, one may ex- 
pect there exists an optimal time Aiop at which the state 
transfer fidelity F gets its maximum value. In Fig. 4 we 
show Xtop versus the intrinsic decoherence rate 7, from 
which one can see that Aiop is shifted to the left-hand side 
of Ato = 7r/2 ~ 1.57, and it decreases with the increase of 
7. Our numerical results also revealed that the magnitudes 
of Atop is independent of the chain length N. 

When considering the average fidelity, the numerical 
results calculated from Eqs. (11), (14), (23), and (24) show 
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tion 2 one can obtain 



Fig. 4. Dependence of Atop on decoherence rate 7. Note that 
the magnitudes of Atop is independent of the chain length A'^. 



" • o<>ooooooooo<l > 




Fig. 5. (Color online) Dynamics of the average fidelity F for 
the XX chain with modulated interactions, where the deco- 
herence rate 7 = 0.15. The black lines from top to bottom 
correspond to chain length A*' = 101, 151, and 201; whereas 
the red lines from top to bottom correspond to chain length 
A'^ = 100, 150, and 200. The inset shows the asymptotic value 
of the average fidelity for infinite time t versus chain length A''. 



that it displays qualitatively the similar behaviors with 
that displayed in Fig. 3. The average fidelity decreases 
with increasing value of both odd and even N , respec- 
tively, and the chain with odd-number qubits seems to be 
more robust on creating high-fidelity state transfer in the 
presence of intrinsic decoherence (see Fig. 5). Moreover, 
as can be seen from the inset of Fig. 5, the average fidelity 
goes to a steady state value in the limit of t — 00, which 
has no relation with the decoherence rate 7. They decrease 
with the increase of both odd and even N, and approach 
the asymptotic value 0.5 in the limit of iV — >■ 00. 

In the absence of intrinsic decoherence (i.e., 7 = 0), 
the above interaction-modulated spin chain can also be 
used to perfectly transfer an entangled state from one 
end of the chain to another [4] . When the decoherence is 
present, however, this ideal spin channel will be destroyed. 
For example. If one start with the Bell state — 
(|01) ± |10))/-\/2 on the first two qubits of the chain, the 
temporal evolution of the concurrence will behaves simi- 
larly as the state transfer fidelity, i.e., it acts as an ampli- 
tude damping channel. When the rescaled evolution time 
\t approaches infinite, from the formulae described in Sec- 



N 



■steady 
N-l.N 



2n 



2n 



(26) 



In fact, one can show that for the initial state \tp) = 
a|01) ± &|10) (|ap + = 1) prepared on the first two 
qubits, the following relation holds 



N 



steady 
N-1,N 



n— 3 



2n- 5 
2n-4' 



(27) 



where Cj"^*'^' = 2|a6| denotes the concurrence of the initial 
state of the first two qubits. This indicates that when the 
rescaled evolution time At approaches infinite, the system 
goes to a steady mirror-symmetric state with two-qubit 
reduced density matrix Pnm(t) = pff-n+i,N-m+i(t), and 
the steady state value C^^'^^i^iN) — Cf'':^'^^{N) decreases 
as the chain length N increases. 

We now investigate entanglement distribution between 
two distant parties through the intrinsic decoherence spin 
channel. For this purpose, we assume the entangled state 
IV') = (|01) -I- 110))/-\/2 is initially prepared between a non- 
interacting qubit NI and the first qubit A on the chain, 
then after some time t, the entanglement will be estab- 
lished between NI and the target spin B. The overall Hamil- 
tonian of the system can be written a.s H' = I (E) H , and 
with the same method used above, one can demonstrate 
that the concurrence Cm,B{N,t) (Note that here N de- 
notes the length of the interacting-spin chain, and does 
not include the noninteracting qubit NI) also behaves as 
a damped oscillation, and when t — ^ 00, we obtain 





(Af+4)/2 



if e odd, 



TT 2n-5 (28) 
II 7 if A^ e even. 



n— 3 



2n-4 



This equation shows clearly that the XX chain with 
even-number qubits is more robust than its counterpart 
with odd-number qubits on distributing quantum entan- 
glement. This is somewhat different from that of the aver- 
age fidelity (see Fig. 5), where the chain with odd-number 
qubits is more efficient on creating high-fidelity state trans- 
fer in the presence of intrinsic decoherence. 



4 Creating entanglement in decoherence 
environments 

In this section, we see intrinsic decoherence effects on 
the creation of entanglement in various kinds of spin net- 
works. For this purpose, we consider the multiarm struc- 
ture M{Ii,12,Na) of the XX Hamiltonian (22) with the 
addition of the exchange couplings between the hub site 
and its nearest-neighbor output sites satisfy the branching 
rule [10]. Here h and I2 denote the number of sites in the 
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Fig. 6. (Color online) Sketch of the multiarm structure of 
the spin network, where the green circle denotes the hub. The 
number of sites in the input and output arms are h = 3, h = 4, 
and the number of output arms is Na = 3. 

input and output arms, respectively, and Na is the num- 
ber of output arms (see Fig. 6). It has been shown that 
in the absence of decoherence environment, this structure 
can be employed to create multi-qubit entangled W state 
at the ends of the outgoing arms. 

The Hamiltonian in the subspace Hi spanned by the 
basis vectors |n) (n = 1, 2, N ) is 

^ = E^«J-I^)01+H-C-, (29) 

(y) 

where the summation runs over all pairs of neighboring 
spins. For the sake of simplicity, we first consider the Y- 
shaped structure Y{li,l2,2). The total number of sites 
now is = /i -I- 2?2 + 1- To examine temporal evolution of 
the concurrence of the prepared initial state, we make the 
following basis transformation for spins just in the same 
position of each arm 

\n^)^±i\n)±\n% (30) 

where n, n' > li + l. Then in the subspace spanned by \n) 
{n ^ li + 1) and |n*) (n > li + 1), the Hamiltonian can 
be rewritten as 

•^n,n+l|")('^+l| + 

n</i + l 

E E Jn..n+i\nn{{n+ir\ + 

n>/i + l ■r=+,- 

V2Ji,+ij,+2\h + l){ih + 2)+| + H.c. (31) 

Clearly, under the transformation (30) the Y-shaped 
structure is transformed into a linear chain consisting of 
the input arm, the hub and one output arm while the other 
output arm is decoupled (see Fig. 7), i.e., this structure 
is identical to the interaction-modulated one-dimensional 
XX chain with chain length I = li + I2 + I- 

For this spin network, if we prepare initial state |1) in 
the first node of the input arm, then after some time t, 
entanglement will be established between the end nodes of 
the output arms (for Fig. 7, it corresponds to node 8 and 
8'). From the formalism described in Section 2 one can 
obtain C{li,l2, X,t) = F{l,X,t). This implies that under 




Fig. 7. (Color online) (a) Sketch of the Y-shaped spin network, 
(b) Under the transformation (30) the Y-shaped structure is 
transformed into a linear chain consisting of the input arm and 
one output arm while the other output arm is decoupled. 

the influence of intrinsic decoherence, the concurrence of 
the created entanglement between the end nodes of the 
output arms also behaves as a damped oscillation as the 
rescaled time Xt evolves. For infinite rescaled evolution 
time Ai — >■ c», the concurrence goes to a steady state value 
C'^^'^yiluh) = 22-2i J]J^^2[4_ 2/(fc - 1)], which can be 
obtained directly from Eq. (25). 

Similarly, for the multiarm structure M{li, I2, Na), us- 
ing the same method, one can obtain that the concurrence 
measuring pairwise entanglement between arbitrary two 
qubits of the end nodes of the output arms is given by 
C{li,l2,NA,X,t) = 2F{l,X,t)/NA (when Na = 2, this 
equality reduces to that describing the Y-shaped struc- 
ture), which observes the similar behaviors as the Y-shaped 
structure, i.e., it behaves as a damped oscillation as the 
rescaled time Xt evolves, and when Xt ^ 00, it goes to a 
steady value C"*°^'^y(/i, ^2, Na) = 2^-^^ nL=2[4 - 2/(fc - 
WNa- 

5 Modified spin chains for high-fidelity state 
transfer 

From the above arguments one can see that the interaction- 
modulated ideal spin channels for perfect state transfer 
are destroyed in the presence of intrinsic decoherence en- 
vironments. Though there exists an optimal rescaled time 
at which one can get a relative high transfer fidelity, how- 
ever, this transfer fidelity (including the average fidelity) 
decreases as the chain length N increases, which puts great 
constrains for long distance communication in interacting- 
spin systems. 

Here we demonstrate that a minor modification of the 
exchange interactions between the first and the last two 
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•^0 \/2x5 ^hx4 \/4x3 ^.hx2 Jo 

(a) o o o — o — o o o 

/o 1 1 1 1 -/o 

(b) o o o o o o o 

Fig. 8. (Color online) Sketches of the two modified seven-site 
spin chains which may serve as spin channels for near perfect 
state transfer in the presence of intrinsic decoherence environ- 
ments. Here we choose A*' odd for it enables more efficient (i.e., 
high speed) state transfer than its counterpart with even A'^. 



nodes of the above structure can fulfill the requirements 
of long distance and near perfect state transfer (see Fig. 
8a). To see this, we display our numerical results for chain 
length = 11 and iV = 51 in Fig. 9(a), from which 
one can see that for all decoherence rate 7, the maxi- 
mum transfer fidelity Fmax approaches unity if Jq is small 
enough (note that when Jq = 0, F = 0), which indicates 
that even in the presence of intrinsic decoherence envi- 
ronments, one can still achieve near perfect transfer of an 
excitation between the opposite ends of a XX chain by 
varying the strength of the exchange interactions between 
the first and the last two nodes of the modulated spin 
chain. 

Another structure which may serve as near perfect spin 
channel for long distance transfer of an excitation in in- 
trinsic decoherence environments is the XX quantum wire 
with the neighboring couplings except those between the 
first and the last two nodes are the same (see Fig. 8b). 
This chain can serve as spin channel for an almost perfect 
state transfer in the absence of decoherence environments 
[11]. When the intrinsic decoherence is present, from Fig. 
9(b) one can see that a long distance transfer of an ex- 
citation whose fidelity can be arbitrarily close to unity is 
also possible for very small but nonzero Jq, even for large 
decoherence rate 7. 

To understand the above phenomenon, we sketch dy- 
namics of Ci{t) in Fig. 10, where Ci{t) = y/ {i\pi{t)\i) de- 
notes the amplitude of the coefficient for the state 
From these two figures one can see that for very small 
but nonzero Jq, except the two end spins (here is node 
1 and 11), the spins at the odd-number sites remain al- 
most unexcited during the time evolution process (i.e., 
Ci(OliGodd;i/i,N — 0), as if the excitation is transferred 
only through the even-number nodes. In fact, a more de- 
tailed analysis show that if Jq — >■ 0+ and t < top, the 
mixedness of the system remain almost unchanged (1 — 
trp"^ ~ 0), i.e., the state here is very close to a pure 
state during the evolution process and thus can be de- 
scribed approximately by \tp{t)) ~ ci{t)\l) + CN{t)\N) + 
I]igevenCi(OI«): whcrc ci{t)+CN{t) ~ 1 (notc that cl{t) + 
c%{t) < 1 since Ci(t)|ig 

even 

In order to better understand why the spins located at 
odd-number sites (the two end spins are exceptions) re- 
main almost unexcited during the time evolution process 
for very small values of Jq , we graph the effects of decreas- 
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Fig. 9. (Color online) Dependence of Fmax on Jo for the two 
modified spin channels displayed in Fig. 8. In both figures the 
black lines from top to bottom correspond to chain length A'^ = 
11 and decoherence rate 7=0.15, 0.30 and 0.45; whereas the 
red lines correspond to A'^ = 51 and 7=0.15. Note that when 
Jo = 0, F = 0. 



ing Jo on the eigenvector population — {\{k\i)\'^}Ek 
for the two modified spin chains in Fig. 11, from which 
one can see that with decreasing values of Jq , the distribu- 
tion of the eigenvectors becomes narrower and narrower. 
Particularly, in the limitation Jq — )■ 0+, {|(fc|l)P}£;fc 
only for the three central eigenvectors of the system, i.e., 
Cfc,i ^ only when ki = {N - l)/2, k2 = (N + l)/2 and 
— {N + 3)/2. Moreover, for these three values of k, it 
can be shown that when Jq — > 0^, Ck.n 7^ only when 
n = 1, N and n € even for k = ki, k^, and n = 1, N for 
k — k2- For example, when = 11, the eigenvectors of 
the Hamiltonian describing the first modified spin chain 
for the above three values of k can be obtained as 

l^i.3>-^(|l> + |iV))Ty5(|2) + |10))T 

yi|6)±y^(|4) + |8)), (32) 

Similarly, for the second modified spin chain, its eigen- 
vectors for the above three values of k with arbitrary chain 
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300 



Fig. 10. (Color online) Dynamics of Ci{t) for the two modified 
spin channels displayed in Fig. 8 with Jo = 0.02, 7 = 0.15 and 
A'^ = 11. Here the curves for even- number i (denoted by blue 
lines) are almost overlapped. 



length N {N G odd) can be obtained as 

f |^o) = ^[|i> + (-i)(^+^'/^|iv)]± 

V-l ^ ^ ' 



2' 



V2| 



(33) 



l^2) = ^[|l) + (-l)(^-^)/^|A^>] 



On the other hand, for initial state |1) prepared in the 
input node A, the density matrix at arbitrary time t can 
be obtained by choosing Q — 11/2, = of Eqs. (10) and 
(11) as 



N 



N 



P{t) = E E Cfc,lCfe'4Cfe^„Cfc',m 
n,m— 1 k,k'—l 



exp 



\n){m\. 



(34) 



Combination of Eq. (34) with the above arguments, 
one can conclude that the spins located at odd-number 
sites remain almost unexcited during the time evolution 
process, i.e., the excitation is transferred only through the 
even-number nodes for very small but nonzero values of 
Jo. 

In fact, from the formalism described in Section 2, one 
can see that the transfer fidelity of an excitation from one 
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Fig. 11. The eigenvector populations {K^|l)l^}Bfc (denoted 
by the relative height of the vertical lines) for the system of 11 
spins. The graphs from left to right in every plot correspond to 
the values of k increases one by one (or equivalently, the eigen- 
values increase one by one). The left seven panels correspond 
to the first modified spin chain with Jo=3, 2.5, 2, 1.5, 1, 0.5 
and Jo — !■ 0^ (from top to bottom), while the right six panels 
correspond to the second modified spin chain with Jo=l, 0.8, 
0.6, 0.4, 0.2 and Jo -> 0+ (from top to bottom). 



end of the chain to another is completely determined by 

O-NN = E '^k,lCk' ,lCk,NCk' ,N 
k,k'—ki 



exp 



^it{Ek-Ek')-^{Ek-Ek'f 



, (35) 



where the coefficients of the eigenvectors for the three val- 
ues of k are given by 



1 

Cfci,l = Ck3,l — -, Cfej^AT — Ck^.N 



(_l)(A'+l)/2 



Ck-i.l — Ck2,N 



(_l)(JV-l)/2 



(36) 



On the other hand, the eigenvalues of the two modified 
spin chains correspond to the three eigenvectors { | A;i ) , | A:2 ) , 
jfcs)} can be written as {—Eo,0,Eq} {Eq can be obtained 
numerically), thus from the above two equations and the 
formalism described in Section 2, one can obtain the trans- 
fer fidelity of one excitation as 

= ^ + J cxp{-2-fE^t) cos{2Eot) - 

i exp (^^^ j cos(i;ot) (TV e odd). (37) 

Since Eq is very small in the limit of Jq — >■ 0^, from 
Eq. (37) one can see that an almost perfect transfer of 
an excitation from one end of the chain to another oc- 
curs after time tc ^ tt/Eq. Moreover, we noted that Eq 
decreases with the decrease of Jq, thus the critical time 
at which the transfer fidelity F gets its maximum value 
increases with the decrease of Jq, and is independent of 
the decoherence rate 7. These conclusions can be corrobo- 
rated by numerical results displayed in Fig. 10. For the sys- 
tem parameters adopted there (i.e., N = 11, Jq — 0.02), 
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Eq can be obtained numerically as £^01 — 0.013801 and 
£■02 ^ 0.012648, thus one has td ~ 227 and td ~ 248. 
Clearly, these results agree well with those displayed in 
Fig. 10. 

A similar analysis shows that for the two modified spin 
chains with even-number qubits, the behavior of the trans- 
fer fidelity is determined only by the two central eigenvec- 
tors of the system in the limit of Jq — t- 0+, i.e., Ck,i 7^ 
only when ki = N/2 and k2 = N/2 + 1. The eigenvectors 
of the Hamiltonian describing both the two modified spin 
chains for the above two values of k can be obtained as 



and 



1 



\k,) = —[\l) + {-lf/'\N)], 
\k) = ^[\l)-{-ir/'\N)], 



(38) 



with the corresponding eigenvalues —Eq and Eq. Combi- 
nation of these with the formalism in Section 2, the trans- 
fer fidelity of one excitation can be obtained as 

F = ^ - ^ exip{-2'yE^t) cos{2Eot) {N e even). (39) 

From Eq. (38) one can see that in the limit of Jo — >■ 0+, 
except the two spins located at the end nodes, all the other 
spins remain almost unexcited during the time evolution 
process, as if the excitation is transferred only between 
the two end nodes. A^Iorcovcr, from Eq. (39) one can sec 
that an almost perfect transfer of an excitation from one 
end of the chain to another occurs after time tc ~ tt/2Eq. 
However, our numerical results show that the values of Eq 
for odd N is much larger than that for even N (e.g., for 
Jo = 0.001, the values of Eq for 11 is about 2524 

and 2804 times larger than that for iV = 10 and 12), thus 
the chain with odd A'' enables a more efficient (i.e., high 
speed) state transfer than its counterpart with even N. 

As a final discussion, it is worthwhile to investigate 
the efficiency of the above two modified spin channels, 
i.e., whether they can serve as near perfect spin channels 
for transfer of an arbitrary one-qubit state by varying the 
strength of Jq. For this purpose, we compute the average 
fidelity. From the above analysis one can obtain straight- 
forwardly that in the limit of Jq 0+, the average fidelity 
can be expressed as 



|6jv| cos(a) (Inn 1 



3 6 2' 

where the coefficients for odd and even A'' are given by 
3 1 



(40) 



ajvjv 



2exp 



(-1) 



exp{-2"fE^t) cos{2Eot) 
~^^ot\ cos(Eot), 

--fE^t 



(A - l)/2 



1 — exp 



cos{Eot) 



Unn = ^~\ exp(-27j;ot) cos(2£;oi), 



bN = (-l)^/2iexp 



-jE^t 



sin{Eot). (42) 



As pointed above, Eq is an infinitesimal in the limit of 

Jo 0+, thus from Eqs. (40) and (41) one can see obvi- 
ously that for odd N, the average fidelity approaches unity 
after time tc ~ tt/Eq. For even N, however, due to the fact 
that a = 7r/2 is not a multiple of 2n, the average fidelity 
can only reach its maximum value 2 /3 (equals to the clas- 
sical average fidelity) after time tc ^ tt/2Eq. This implies 
that for the above two modified spin channels with even 
A'', one cannot achieve near perfect state transfer of an 
arbitrary one-qubit state simply by varying the strength 
of Jq. But if one can apply an external magnetic field B 
along the z axis of every spin (this does not change the 
eigenvectors of the system since [H, a^^^] = 0), the phases 
of the received state at the destination node may be cor- 
rected. With this method, we performed numerous calcu- 
lations and the numerical results revealed that the average 
fidelity can also approaches unity by choosing appropriate 
strength of the magnetic field (e.g., for A^ = 10, Jo = 0.01 
and B = ±0.0003, the two modified spin channels give 
rise to F^ax = 0.9998 and F^ax = 0.9996, respectively). 



6 Summary 

To summarize, we have investigated quantum state trans- 
fer, generation and distribution of entanglement in the 
model of Milburn's intrinsic decoherence environment. We 
focused on diverse interaction-modulated spin networks 
which may serve as perfect spin channels in the absence of 
decoherence. As one expected, the state transfer fidelity 
as well as the amount of the generated and distributed 
entanglement will be significantly lowered by the intrinsic 
decoherence environment, and this detrimental effects be- 
come severe as the decoherence rate 7 and the spin chain 
length A^ increase. For infinite evolution time t, we show 
analytically that both the state transfer fidelity (including 
the average fidelity) and the concurrence of the generated 
and distributed entanglement approach steady state val- 
ues, which are independent of the decoherence rate 7. This 
brings great constraints on these structures as spin chan- 
nels for long distance and high-fidelity communication. Fi- 
nally, as alternative schemes to diminish the detrimental 
effects, we presented two modified spin chains which may 
serve as spin channels for long distance and near perfect 
state transfer in the intrinsic decoherence environments. 
Our results revealed that in the limit of Jq — >■ 0+, these 
two modified spin channels generate maximum fidelity 1 
after certain time tc ~ tt/Eq for spin chains with odd- 
number qubits. For spin chains with even-number qubits, 
however, one needs to apply an external magnetic field in 
order to achieve near perfect state transfer. 
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